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Abstract: In this study, a new method for the estimation of aeolian sand transport rate was 
developed; the method employs a ceramic sand flux sensor (UD-101). UD-101 detects 
wind-blown sand impacting on its surface. The method was devised by considering the 
results of wind tunnel experiments that were performed using a vertical sediment trap and 
the UD-101. Field measurements to evaluate the estimation accuracy during the prevalence 
of unsteady winds were performed on a flat backshore. The results showed that aeolian 
sand transport rates estimated using the developed method were of the same order as those 
estimated using the existing method for high transport rates, i.e., for transport rates greater 
than 0.01 kgm"^ s~\ 

Keywords: piezoelectric sensor; impact counts; wind-blown sand transport; instantaneous 
sand transport; moisture content 



1. Introduction 

Aeolian sand transport is an important process that is studied in many flelds such as physics, 
geosciences, agriculture, and coastal engineering. Considerable effort has been devoted to 
understanding sand transport mechanisms, including the aerodynamic entrainment of sand grains, 
acceleration of sand grains in air streams, and transport of sand by steady winds [1-16]. 

Aeolian sand transport mechanisms observed in the field under unsteady wind conditions differ 
signiflcantly from those observed under steady wind conditions. Since the 1990s, high-frequency 
sampling systems such as piezoelectric mass flux sensors [17-22], high-temporal-resolution sediment 
traps [23-27], and microphonic flux sensors [28] have been employed for the detection of aeolian sand 
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flux at frequencies of several hundreds of hertz and higher. Wind tunnel and field experiments 
performed using these systems have improved our understanding of sand transport mechanisms under 
unsteady wind conditions. However, further investigations are required, especially under complex 
field conditions. 

The ceramic sand flux sensor UD-101 (Chuo Kosoku), which is a piezoelectric sand flux sensor that 
can measure the flux of aeolian sand impacting on its surface, was developed by Kubota et aL [29,30] 
(see Figure 5c). The sensor detects piezoelectric signals generated by collisions between aeolian sand 
grains and the sensor surface using the same principle as that employed in Sensit [17] and Safires [19]. 
A significant difference between UD-101 and the other piezoelectric sensors is that UD-101 is 
unidirectional and does not self-orient toward the wind direction, while both Sensit and Safires are 
omnidirectional. However, UD-101 has the following advantages over the other piezoelectric sensors: 

(i) The sensor is calibrated by performing both wind tunnel and field experiments. 

(ii) The sensor records an impact value of 0 or 1 at 10 kHz; it samples once per second, and 
theoretically, it can detect a maximum of 10,000 impact counts in the sampling interval, 
i.e., the sensor has high temporal resolution. 

(iii) The sensor is relatively inexpensive. 

Udo et al. [22] carried out short-term field observations at the Hasaki Beach in Japan using UD-101 
during the period 12-16 January 2005. The following features were inferred from the UD-101 data: 

(i) During periods without rainfall and in the presence of longshore winds (conditions 
similar to those employed in wind tunnel experiments), the aeolian sand fiux at a certain 
height above a fiat ground surface increased significantly with the wind velocity 
and approximately equaled the fiux estimated in wind tunnel experiments using 
Kawamura's [3] equation. 

(ii) The fiux decreased significantly with an increase in precipitation, i.e., with an increase in 
the moisture content of the sand surface; however, even during periods with rainfall, fiux 
was detected during strong wind conditions. 

(iii) The fiux increased with a decrease in the angle between the sensor direction and 
wind direction. 

Further, Udo [31] carried out relatively long-term observations during the period from April 2005 
to January 2006 at the Hasaki Beach and showed that a large amount of sand was transported during 
typhoons. This result implies that UD-101 can be employed for performing long-term fiux 
measurements even during storm events. As the next step toward making long-term flux measurements 
in fields possible, it is necessary to develop a measurement method for not only the aeolian sand flux 
at a certain height above the ground surface but also the aeolian sand transport rate (e.g., the total 
aeolian sand flux between the ground surface and the maximum height of the aeolian sand layer). In 
particular, the accurate estimation of the rate of seasonal, yearly, or decadal sediment transport due to 
both littoral drift and aeolian transport at beaches would be invaluable for coastal management. 
However, the aeolian transport rates estimated by most researchers have not been validated because of 
the lack of practical estimation/measurement methods; examples of such estimations are estimations 
made by using conventional equations along with the mean wind data [32] or by integrating areas of 
beach morphological changes from a boundary under the assumption that the sand transport rate at the 
boundary is zero [33]. 
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The objective of this study was to develop a new method employing UD-101 for the estimation of 
the aeolian sand transport rate [29,30]. We applied the method to field data obtained during the period 
from 22 September to 6 December 2005 [31] and assessed the estimation accuracy to examine the 
suitability of the method for the direct estimation of long-term transport rates in the field. 

2. New Estimation Method for Determining Total Aeolian Sand Flux 

2.1. Characteristics of UD-101 Data 

Figure 1 shows a time series of the instantaneous horizontal wind velocity u and blown-sand 
impact count n measured at the Hasaki Beach; the median sand grain size Jso at the beach 
was 0.2 X 10 m [22,31]. n shows a positive relationship with u during periods without any rainfall, 
but decreases significantly after a rainfall event. The sensor cannot detect the impact of blown sand if 
wet sand adheres to its surface after a rainfall event; however, it can detect the impact if this sand is 
subsequently removed by heavy rainfall. The sensor does not detect the impact of raindrops. 

Figure 1. Time series of instantaneous horizontal wind velocity u and blown- sand impact 
count n at the Hasaki Beach [22,31]; the data were recorded for 5 min at a frequency 
of 1 Hz (a) during intervals without rainfall on 12 January 2005 and (b) after a rainfall 
event on 25 September 2005. The wind velocity was measured using an ultrasonic 
anemometer at a height of 0.9 m above the fiat ground surface {zu = 0.9 m). The impact 
count was measured using UD-101 at = 0.04 m. The horizontal distance between the 
wind and saltation sensors was 1 m. The median sand grain size dso of the ground surface 
at the beach was 0.2 x 10 m. 
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Many studies have shown that the vertical wind velocity profile in a wind tunnel is logarithmic [34,35]. 
In previous studies, the threshold wind velocity Ut has been estimated from both a wind profile law and 
Bagnold's [1] equation. The logarithmic wind profile law is expressed as: 

M = — ln(z„/zoJ (1) 

Here, m is the wind shear velocity, /ris the von Karman constant (= 0.4), Zu is the height at which u 
is measured, and z^s is the saltation roughness length. Saltation is the hopping motion of sand grains 
and the primary mode of aeolian sand transport in beach intertidal and supratidal zones. Aeolian sand 
transport leads to an increase in the aerodynamic roughness length. Ut is obtained by substituting 
m = u*t in Equation (1). The threshold wind shear velocity is estimated from 

Bagnold's [1] equation: 



u.,=a^^-^gd (2) 

where a is the numerical transport rate coefficient (= 0.1), is the air density, cris the sand density 
(= 2,650 kg m~^), g is the gravitational acceleration, and d is the sand grain size, zo^ is estimated from 
the relation proposed by Chamock [36] and Chamberlain [37]: 

2 

^0. =^0— (3) 
where co is a constant (=0.16 for saltation on beaches [38]). 

2.2. Relationship between Impact Counts and Aeolian Sand Flux at a Measurement Height 

Kubota et aL [30] performed experiments in a 20-m-long wind tunnel using the experimental setup 
that was employed by Hotta et al. [39]; the sand used in the experiments had a d^o value 
of 0.25 X 10 m. The authors showed that the impact counts approached a constant value 
asymptotically at frequencies greater than 8 kHz. They also compared the sand mass flux qs measured 
in the wind tunnel using an 8-kHz sensor [30] with q measured using a vertical-distribution-type sand 
trap [39]. qs at the measurement height Zs was estimated from n using the following equation [22]: 

/ \ lod.J'n 

<ls[^s)—rTT- (4) 

where ^ is the sampling period oi n {= \ s) and ds is the diameter of the sensor (= 0.012 m). The 
vertical-distribution-type sand trap was developed by Hotta and Horikawa [40], and it can be used to 
measure the flux at 37 different heights in the range 0-0.505 m. q was corrected for trap efficiency. 
The wind tunnel experiments performed by Kubota et aL [30] and Hotta et aL [39] showed that qs{zs) 
for (iso = 0.25 x 10"^ m was related to q{zs) measured using the vertical-distribution-type sand 
trap [22], which was closely approximated as: 

q{z^) = \Q.51qXz,) (5) 
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2.3. Relationship between Aeolian Sand Flux at a Measurement Height and Total Sand Flux 

The relationship between the sand flux at z^, q{zs) [kg s~^], and the total sand flux Q [kg s~^] 
is shown in Figure 2. Hotta et al. [39] employed the same wind tunnel as that used by 
Kubota et al. [30], with 0.17 m s"^ < m* < 1.77 m s"' and 0.14 x IQ-^ m < J50 < 0.68 x IQ-^ m. Here, 
only the data for (iso = 0.25 x 10 m are shown. Q was found to be approximately proportional to q{zs), 
as expressed by: 

Q = M{zs) (6) 

where the regression coefficient R of the approximation was 0.96, 0.99, 1.00, 1.00, 0.99, 0.98, 0.97, 
and 0.97 for z, = 0.005, 0.02, 0.04, 0.06, 0.08, 0.10, 0.12, and 0.14 m, respectively. As shown in 
Figure 3, the constant A shows a positive relationship with z^, and it can be expressed as: 

^ = 5exp(CzJ (7) 

The relationships between z^ and ^ and between Zs and A (shown in Figures 2 and 3, respectively) 
indicate that Zs should be set in the range 0.02-0.08. 

Figure 2. Relationship between the sand flux at z^, q{zs), and the total sand flux Q for 
^/50 = 0.25 X 10"^ m. 
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Figure 4. Relationship between (a) dsQ and B and (b) dsQ and C. 
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The constants B and C show a positive relationship with Jso (Figure 4): 



and: 



5 = 89.7^50+0.0317 



C = -1.21xlOV5o+20.4 



Finally, the following equation is derived from Equations (4)-(9): 



0.0006 



Q = 



- (89.7^/50 + 0.03 1 7)exp[(- 1 .2 1 X 1 0' d,^ + 20a)z^ ] 



0.0008 



(8) 
(9) 

(10) 



This equation suggests that Q can be estimated only if n, dso, and Zs are measured by performing 
field experiments, although the equation was derived from observations made during 20-m-long wind 
tunnel experiments performed by Hotta et al. [39] and Kubota et al. [30] on a flat ground surface under 
steady wind conditions with < 1.8 m s~^ and at a sensor height of 0.005 <Zs< 0.14. 



3. UD-101 Data Obtained at Hasaki Beach 



3.1. Summary of Field Measurements 

Field measurements were performed at the Hasaki Beach using UD-101. The beach faces the 
Pacific Ocean and is located in Ibaraki Prefecture [35°50'25''N, 140°45'42''E; Figure 5(a)] [22,31]. 
The beach is of the dissipative type and has a foreshore slope ranging from 0.02 to 0.05 and a width 
of 100 m [Figure 5(b)]. ds{) is approximately 0.2 x 10"^ m along the backshore. Beach grasses grow 
around a small embryo dune at j; = 50 m and around a longshore dune at j; = 70 to 90 m just in front of 
the coastal forest area, where y is the onshore distance from the average shoreline position in 2004. 
The measurement point was located in a flat area between the embryo dune and backshore dune. 

The annual wind direction is predominantly in the north-northeast direction [Figure 5(a)]. Winds 
are relatively weak in summer, except during typhoons, but strong in winter due to the presence of 
low-pressure systems, u, the wind direction 6, and n [the same data obtained by Udo [31], Figure 5(c)] 
were measured from 21 September to 3 December 2005. u and ^were measured at a height of 0.9 m 
above the measurement point (z = 0.9 m) using a three-axis ultrasonic anemometer (Delta Ohm; 
model: HD2003), while n was measured at z = 0.04 m using UD-101. Four UD-101 sensors were 
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deployed, and they were directed toward the north, east, south, and west. The wind direction in the 
horizontal plane is expressed clockwise from north, from 0° to 360"^. The data of u, 0, and n were 
recorded on a data logger system (FieldPoint cFP-2000, National Instruments) at a frequency of 1 Hz 
and were automatically downloaded on a weekly basis from the logger. The precipitation Pr was 
measured at 10-min intervals at the Choshi Meteorological Observatory (35''44'18''N, 140''5r24''E) by 
the Japan Meteorological Agency. 

Figure 5. (a) Location of the study site (Hasaki Oceanographical Research Station, HORS) 
and the wind rose for 2004. (b) Topography of the Hasaki Beach and the location of the 
measurement point. The cross-shore distance to the average annual shoreline is 0 m; the 
offshore direction was considered positive, (c) Location of the measurement point 
on 17 January 2005. The ultrasonic anemometer and ceramic sand flux sensors were 
installed as shown. The four UD-101 sensors were oriented toward the north, east, south, 
and west directions. 





In addition, in this study, Zs was measured manually almost every day in order to calculate Q using 
Equation (10). The instantaneous Zs was calculated by linear interpolation. The Zs value of the sensor 
oriented toward the wind direction tended to be large and that of the sensor directed in the opposite 
direction tended to be small; for example, for north winds, the Zs value of the sensor directed 
northward tended to be large and that of the sensor directed southward tended to be small. 

The measurement results are shown in Figure 6. Data could not be obtained from 28 October 
to 13 November 2005 because of network problems. Large amounts of sand were transported during 
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the transport events on 24-26 September due to the impact of typhoon 0517, which resulted in a total 
precipitation of 39.5 mm; similar events involving the transport of a large amount of sand occurred 
on 11-12 October (period without rainfall); on 18-20 October (due to the impact of typhoon 0520), 
when the total precipitation was 6 mm; and on 14-15 November (period without rainfall) [31]. 

Figure 6. Time series of 10-min averages of the horizontal wind velocity Umeam horizontal 
wind direction Omeam blown-sand impact count rimeam estimated total aeolian sand flux 
Qmean, aud prccipitatiou Pr; the time series of Qo estimated from Equation (13) and the 
daily sensor height from the ground surface Zs are also shown. The solid gray lines in the 
top figure show the maximum and minimum instantaneous wind velocities. The wind 
direction is expressed clockwise from north, from 0"^ to 360°. Only the maximum rimean for 
the four directions in each of the 10-min intervals is plotted. 
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3.2. Evaluation of the Accuracy of the New Method by Comparing Q with Qb, Qk, Qo, cind Ql 



The Q value estimated from n was compared with that obtained using the conventional equations 
for a constant wind velocity in order to examine the accuracy of Q estimated from the UD-101 data. 
Here, the commonly employed equations of Bagnold [1], Kawamura [3], Owen [5], and Lettau and 
Lettau [6] (for Qb, Qk, Qo, and Ql, respectively) were used [13,16,41-43]: 
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(11) 



(12) 



Qo 



Cq — M* (m* - + ) 



(13) 



2i =^^M— J (14) 
VZ) g 

Here, D is the reference sand grain diameter (= 0.25 x 10 m) and cb, ck, cq, and cl are constants 




with values of 1.8, 2.78, 0.25 + wo/3i^* (where wo is the fall velocity), and 4.2, respectively. The sand 
transport rate was estimated to be zero when < was estimated from Equation (1) using the 

measured value of w; was obtained using Equation (2) and wo was determined from 
Rubey's [44] equation: 



where s is the relative density of the sand (= dp) and v is the kinematic viscosity of air. These 
equations were derived for a constant wind velocity, and hence, Q estimated using the UD-101 data 
was averaged over 10 min. The 10-min mean sand flux Qmean during periods without rainfall was 
compared to Qb, Qk, Qo, and Ql. 

Figure 7 shows a comparison of Qmean values obtained from n with Qb, Qk, Qo, and Ql values 
derived from the equations given above, for periods without rainfall. Data for w = 0 and m < un 
(i.e., u < 4.05 m s~^) were excluded. Qmean showed a clear positive relationship with Qb, Qk, Qo, and 
Ql, and the estimated values followed the order Qo< Qb< Ql< Qk- The best agreement was observed 
between Qmean and Qo when 0.02 < Qmean < 0.10 and between Qmean and Qb when Qmean > 0.10. 

Figure 8 shows double logarithmic plots of Qmean versus Qb, Qk, Qo, and Ql for east-northeast 

(15° < 0mean < 105°), SOUth-SOUthcast (105° < 0mean < 195°), WCSt-SOUthwCSt (195° < 0mean < 285°), and 

north-northwest (285° < 0mean < 375°) winds. For small Qmean values, i.e., less than 0.01 kg m~^ s~\ the 
Qmean valuc was Overestimated by Qb, Qk, Qo, and Ql. The probable causes for the overestimation are 
the low measurement accuracy of n for small Qmean, the low estimation accuracy of zos for small u*, 
variations in the wind velocity, and the high moisture content at the beach surface. Raupach [38] 
showed that cq in Equation (3) depended on u*, especially for small w*. A decrease in the instantaneous 
wind velocity below the threshold velocity, even when the mean wind velocity is greater than the 
threshold velocity, can lead to an overestimation. The moisture content of the backshore is nonzero 
even during periods without rainfall because of the proximity of the backshore to the sea. 
Udo et al. [22] demonstrated that the sand transport flux at a height, qmean, decreased by more than one 
order of magnitude when the moisture content increased during rainfall periods. 




(15) 
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Table 1 lists the linear regression relationships between Qmean estimated from n and Qb, Qk, Qo, and 
Ql for the east-northeast, south- southeast, west-southwest, and north-northwest winds (see Figure 8). 
While a strong correlation was observed for the north-northwest wind (in the longshore direction), a 
clear correlation was not observed in the case of the other winds because only a small number of Qmean 
values were greater than 0.01 kg m~^ s~\ Udo [31] showed the dependency of qmean on dmean and that 
gmean was high whcu dmean was iu the lougshorc direction (especially in winter); however, qmean was 
low when dmean was seaward (in late summer and winter) and was even lower when dmean was 
landward (in spring and autumn). The differences between Qmean and Qb, Qk, Qo, and Ql were found 
to be independent of dmean- 



Figure 7. Comparison of Qmean estimated from n with Qb, Qk, Qo, and Ql estimated using 
the equations given in the text, when ni^O or m> u*t (i-e., u > 4.05 m s~^) during periods 
without rainfall. The solid line is the linear regression line. The dashed line is the line that 
would be obtained in case of perfect agreement. 
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Figure 8. Comparison of Qmean estimated from n with Qb, Qk, Qo, and Ql for 
east-northeast (15"^ < Omean < 105"^), south- southeast (105° < Omean ^ 195°), west-southwest 
(195° < Omean < 285°), and north-northwest (285° < Omean < 375°) winds. Qb is either zero 

O 1 1 

or greater than 1.76 x 10" kg m s because Equation (11) does not include the term 
(t/* - u^t)- The dashed line denotes perfect agreement. 
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Table 1. Linear regression expressions for the relationships between Qmean estimated from 
n, and Qb, Qk, Qo, and Ql for east-northeast, south-southeast, west-southwest, and 
north-northwest winds. is the regression coefficient and is the number of data points. 
See Figure 8 for related plots. 
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Hotta and Horikawa [40] performed field experiments at beaches by using a trench-type trap and a 
vertical-distribution-type trap to estimate the aeolian sand transport rate and showed that the possible 
ranges of cb and ck are 0.5-2.0 and 1.0-3.0, respectively. This result suggests that the actual values of 
Qb and Qk are possibly one-third of those shown in Figures 7 and 8. All the field results showed that 
the sand transport rate estimated using the new method together with the UD-101 data was of the same 
order as that estimated using the conventional equations, although several differences were 
observed between the aeolian sand transport rate in the wind tunnel and that in field 
experiments [38,40,45,46]. The new method is useful for the direct estimation of the sand transport 
rate in a field analysis. Further analysis using the UD-101 data can provide additional insights into 
aeolian sand transport mechanisms in fields. 

4. Conclusions 

In this study, we developed a new method for the estimation of the aeolian sand transport rate and 
demonstrated the usefulness of UD-101 for field measurements of the aeolian sand transport rate. The 
transport rates estimated from the UD-101 data were reasonable and showed a clear positive 
relationship with the rates estimated using the conventional equations. Further, for 
Qmean ^ 0.01 kg m"^ s"\ the transport rates estimated using the UD-101 data were of the same order as 
those estimated using the conventional equations. The best agreement was observed between Qmean and 
Qo when 0.02 < Qmean < 0.10 and between Qmean and Qb when Qmean > 0.10. It is necessary to develop 
a method for accurately determining the constant values in the conventional equations. The cause for 
the disagreement between Qmean and Qb, Qk, Qo, and Ql for small Qmean should be investigated in a 
future research. In field analyses, the estimation of large transport rates under strong winds is 
important and the proposed estimation method is expected to be useful in obtaining the transport rate. 
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